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In this paper we explore the potential multistability of the climate for a planet around the habitable zone. We focus on 
conditions reminiscent to those of the Earth system, but our investigation has more general relevance and aims at presenting 
a general methodology for dealing with exoplanets. We describe a formalism able to provide a thorough analysis of the 
non-equilibrium thermodynamical properties of the climate system and explore, using a a flexible climate model, how such 
properties depend on the energy input of the parent star, on the infrared atmospheric opacity, and on the rotation rate of the 
planet. We first show that it is possible to reproduce the multi- stability properties reminiscent of the paleoclimatologically 
relevant snowball (SB) - warm (W) conditions. We then characterise the thermodynamics of the simulated W and SB 
states, clarifying the central role of the hydrological cycle in shaping the irreversibility and the efficiency of the W states, 
and emphasizing the extreme diversity of the SB states, where dry conditions are realized. Thermodynamics provides the 
clue for studying the tipping points of the system and leads us to constructing empirical parametrizations allowing for 
expressing the main thermodynamic properties as functions of the emission temperature of the planet only. Such empirical 
functions are shown to be rather robust with respect to changing the rotation rate of the planet from the current terrestrial 
one to half of it. Furthermore, we explore the dynamical range of slowy rotating and phase locked planet, where the length 
of the day and the length of the year are comparable. We clearly find that there is a critical rotation rate below which 
the multi- stability properties are lost, and the ice-albedo feedback responsible for the presence of SB and W conditions 
is damped. The bifurcation graph of the system suggests the presence of a phase transition in the planetary system. Such 
critical rotation rate corresponds roughly to the phase lock 2:1 condition. Therefore, if an Earth-like planet is 1:1 phase 
locked with respect to the parent star, only one climatic state would be compatible with a given set of astronomical and 
astrophysical parameters. These results have relevance for the general theory of planetary circulation and for the definition 
of necessary and sufficient conditions for habitability. 
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1 Introduction 

The investigations of extra- solar planetary objects is one of 
the most active fields of research in astrophysical sciences. 
After the first discoveries dating back to the mid 1990s, im- 
proved observational datasets and techniques of data analy- 
sis have made possible to catalogue and studies hundreds of 
planets orbiting parent stars other than the Sun. These plan- 
ets come in a great variety of physical, chemical, and orbital 
properties. Factors of great relevance include the presence 
or lack of a rocky core, the composition of the atmosphere, 
the color intensity of the radiation emitted by the parent star, 
the orbital parameters, the presence or lack of a tidal-lock 
condition. An especially important problem for rocky plan- 
ets (the so-called super-Earths) is determining under which 
conditions one can expect that the planet may feature the 
presence of liquid water at surface, as this is seen as neces- 
sary condition for the presence of life. We refer the reader 
to some recently published books (Dvorak 2008; Kastings 
2009; Seager 2010; Perryman 2011) . 

The existence of great variety of planetary conditions 
is, on one side, obviously a scientific gold mine, and on the 
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other side it is a maze one can easily get lost in. It is incredi- 
bly attractive to have the possibility of studying the individ- 
ual properties of a fast growing number of newly-discovered 
astronomical objects, for which a growing amount of infor- 
mation is coming and is very likely to come in the near and 
medium-term future, considering the technological devel- 
opment and the coming astrophysical big science initiatives. 
On the other side, one may demand what is in the long term 
the scientific merit of performing a Linnaean compilation 
of the properties of distant planets, and of developing ra- 
diative models and/or general circulation models aimed at 
describing the radiative and dynamical properties of their 
atmospheres. We know that, even for planets we know well, 
like those of our solar system, developing accurate numer- 
ical models able to account for the dynamics and thermo- 
dynamics of the planetary atmosphere compatibly with our 
observations is extremely challenging. And, coming back 
to Earth, it is well-known that there is no satisfactory the- 
ory of climate dynamics, able to account for the complex, 
multi- scale interplay between radiative forcings, dynami- 
cal instabilities, positive and negative feedbacks, impacts of 
slow modulations of orbital and geological conditions. The 
development of a complete climate theory is, in fact, one of 
the great challenges of contemporary science, and the devel- 
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opment of seamless models able to account for the behavior 
of the fluid component of our own planet over a large range 
of temporal and spatial scale is still a distant hope. The the- 
oretical interpretation of the onset and decay of many well- 
documented events of our planet's history - like the ice ages 
- is far from settled, let alone the presence of models able 
to simulate them satisfactorily. A provocative question one 
may formulate using a reductionistic point of view may well 
be: what is the hope/goal/scope of studying the properties of 
the atmospheres of distant planets if we still know relatively 
little of the properties of our own? 

We argue here that, in fact, the investigation of exo- 
planets creates an outstanding scientific horizon for trans- 
disciplinary research across astrophysics and geophysics, 
and, in fact, geophysical sciences can definitely contribute 
to the investigation of non-terrestrial objects, and, on the 
other side, by studying other planets, we can understand 
better the planet Earth, because we can take a broader sci- 
entific perspective, less centered on trying to predict and ex- 
plain the variability of some specific physical properties of 
the Earth's fluids (typically of direct interest for mankind) 
preferentially over other ones. The challenge is then try- 
ing to develop general physical and mathematical theoret- 
ical methods and diagnostic tools to be used for studying a 
large variety of planetary atmosphere, and, at a more practi- 
cal level, to construct flexible numerical models able to sim- 
ulate many different planetary atmospheres once some basic 
input parameters are given. In this perspective, one must re- 



fer to recent paper by |Read (201 1)| where, along the lines of 
what traditionally done in fluid mechanics, it is suggested to 
construct several dimensionless numbers and use them for 
classifying the circulation regimes of classes of planets. 

In this paper we try to give a limited, temptative but 
three-fold contribution to establishing such a link between 
geophysical and astrophysical sciences. We have drawn 
from some recently published material (Lucarini 2009; Lu- 
carini, Fraedrich & Lunkeit 2010a; 2010b; Boschi et al. 
2013) and we add some new results of distinct astronomical 
flavor. First, we introduce some ideas and methods derived 
from basic results of macroscopic thermodynamics that al- 
low to define in a rigorous way indicators describing the 
basic non-equilibrium properties of a general planetary at- 
mosphere: its ability to transform available potential energy 
into kinetic energy, thus performing work like a thermal en- 
gine; and its production of entropy through irreversible pro- 
cess; and how these two concepts can be related in a com- 
pact conceptual framework. Then, we use the derived ther- 
modynamical indicators to provide a novel description of 
a classical problem of paleoclimatology, i.e. the onset and 
decay of the so-called Snowball Earth. Based on the evi- 
dence supported by Hoffman et al. (1998) and Hoffman and 
Schrag (2002)1 it is expected that the Earth is potentially 
capable of supporting multiple steady states for the same 
values of some parameters such as the solar constant and 
the concentration of carbon dioxide, which directly affect 
the radiative forcing. Such states are the presently observed 



state, and a state where the planet is entirely ice covered 
and the surface temperature is much lower than the present 
one. Such conditions hardly allow for the presence of life, 
so this issue is of extreme relevance of the general quest 
for defining habitability condition in other planets. Initial 
research using simple 0-D models (Budyko 1969; Sellers 
1969), 1-D models (Ghil 1976) as well as more recent anal- 
yses performed using complex 3-D general circulation mod- 
els (Marotzke & Botztet 2007; Voigt & Marotzke 2011; 
Pierrehumbert et al. 2011), provide support for the exis- 
tence of such bistability for certain range of parameters of 
the system. The main mechanism triggering the abrupt tran- 
sitions between the snowball and the warm state is the posi- 
tive ice-albedo feedback (Budyko 1968; Sellers 1969). Such 
a feedback is associated with the fact that as temperatures 
increase, the extent of snow and ice cover decreases thus 
reducing the albedo and therefore increasing the amount of 
sunlight absorbed by the Earth system. Conversely, a neg- 
ative fluctuation in the temperature leads to an increase in 
the albedo therefore reinforcing the cooling. We show how 
thermodynamics allows for a much deeper physical inter- 
pretation of the ice-albedo feedback, which is the instabil- 
ity mechanism mainly responsible for the multi stability of 
the system. At this regard, we use the numerical evidence 
gathered by running PlaSim, a general circulation model of 



intermediate complexity (Fraedrich et al. 2005) and study 
some basic properties of the climate states realized when the 
solar constant is modulated between 1160 Wm -2 and 1510 
Wm -2 and the values of [CO2] are varied between 90 and 
2880 ppm, and define the region of multistability. Then, we 
extend the previous result in the direction of providing in- 
formation useful for general planetary objects by altering a 
critical astronomical parameter, i.e. the rotation rate of the 
planet, and explore the case of slowly rotating planets. We 
perform these investigation with a modified climate model, 
where, for sake of generality, the land-sea contrast typical 
of our planet is removed and an ocean only surface is con- 
sidered. Whereas for the climate of a so-called Aquaplanet 
is indeed different from that of the corresponding realistic 
Earth, its structural properties in terms of multi stability are 
the same (Voigt e t al. 2011)1 so we gain generality without 
losing comparability of the results. We then study how the 
structural properties of the system change when we change 
the rotation rate of the planet. We discover that when the ro- 
tation rate become smaller than a critical but non-zero value, 
the multi stability condition is lost, and, in particular, the 1 : 1 
phase locked planet has a unique domain of attraction for its 
physical conditions. The exploration of very diverse dynam- 
ical regimes determined by the rotation rate is performed 
using a prototypical version of a flexible modeling suite 
which is being built starting from the (already rather flex- 
ible) PLASIM platform. Such modeling suite, once com- 
pleted, will allow for simulating extremely different super- 
Earths, using exactly the same dynamical core and modeling 
platform. 
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The paper is organized as follows. In Section [2] we pro- 
vide a basic introduction to some tools of non-equilibrium 
thermodynamics used throughout the paper. In Section [3] 
we describe our numerical simulations and the set-up of the 
model used in this study. In Section [4] we discuss the bista- 
bility properties of the climate state, discussing the Snow- 
ball and the Snow-free states, and the transitions between 
the two. in Section [5] we present the results obtained on the 
role of the rotation rate in defining the structural properties 
of the system. In Section [6] we draw our conlusions. 

2 Non-equilibrium Thermodynamics of the 
climate 

In this section we briefly introduce the thermodynamic diag- 
nostic tools used in this review to analyse the multistability 
properties of an Earth-like planet. An extended discussion 
of such a formalism and of its applications to study global 
properties of non-equilibrium systems can be found else- 
where (Fraedrich & Lunkeit 2008; Lucarini 2009; Lucarini 
et al. 2010; 2010b). 

The total energy budget of the atmosphere can be writ- 
ten as E = P + K, where K = j v pv 2 /2dV represents 
the total kinetic energy and P = j v (c p T + gz + Lq)pdV 
is the moist static potential energy (which is the sum of sen- 
sible, potential and latent energy) and V is the atmospheric 
domain. It can be shown (Peixoto & Oort 1992) that: 



K = —D + W, 
P = * + D - W, 



(1) 

(2) 



where D = J v dVpe 2 is the dissipation of kinetic energy 
due to friction (e 2 > is the local rate of heating associated 
with kinetic energy dissipation), W = -JyV- VppdV 
is the instantaneous work done by the system and ^ = 
f v dV{— V • H) is the total heating due to convergence 
of sensible heat, latent heat and radiative fluxes. Equations 
{TJ-H imply that E = 4/ and therefore the frictional heat- 
ing e 2 does not increase the total energy since it is just an 
internal conversion between kinetic and potential energy. 
Taking the climate as a non-equilibrium steady state sys- 



tem (NESS, see Galla votti 2006| ), we have that over long 

time scales E = P = K = (here and in the follow- 
ing the bar indicates averaging over long time periods) and 

therefore ^ = 0. If we define the total diabatic heating 
Q = pe 2 — V H and split the domain V into the subdomain 
V+ , where Q = Q+ > 0, and V~, where Q = Q~ < 0, it 
can be easily seen from Eq. ^ that: 



^ + L> = P + W 



/ dVpQ + + [ dVpQ- 
Jv+ Jv- 



(3) 

with i>+ > and i>~ < by definition. From eq. |l[ 
U and (|]) it is straightforward to show that D = W — 



4>+ + i> - > 0, which summarizes the Lorenz (1967) energy 
cycle . Therefore the atmosphere can be interpreted as a heat 
engine, in which 4> + and <&~ are the net heat gain and loss 
needed in order to produce mechanical work W given by 
their (algebraic) sum. The efficiency of the atmospheric heat 
engine, i.e. the capability of generating mechanical work 
given a certain heat input, can therefore be defined as: 



rj = ($'+ + $-)/($+) = W/&+. 



(4) 



The analogy between the atmosphere and a (Carnot) heat 
engine can be pushed further if we introduce the total rate 
of entropy change of the system, S = j v Q/TpdV and use 
the V ± partition of the domain: 



S : 



v+ 



pdV 



pdV 



(5) 



where S + > and E < at any time. In a steady state, 
!!+ = — as 5 = and the following expression holds 



6+ 
where 



(6) 



3>± = t ± G ± (7) 

and the temperatures 6 + and 6~ are the time and space 
averaged temperatures of the V + and V~ domains respec- 
tively. Since = |E"| and |$ + | > |$~|, it can be 
shown that 9 + > 0~, i.e absorption typically occurs at 
higher temperature than release of heat (Peixoto & Oort 
1992; Johnson 1997). From Eqs. gj, ^ and we de- 
rive that r] = (6+ - 6")/6 + (Johnson 2000; Lucarini 
2009), similarly to the definition of the actual Carnot effi- 
ciency (Fermi 1956) . 

In a planet the non-equilibrium steady state is main- 
tained by the global compensation between the net radiative 
heating in the warm regions (low latitudes on Earth) and net 
cooling in the cold regions (high latitudes on Earth). The en- 
ergy budget of the warm and cold regions is closed thanks 
to the presence of large scale transports performed by the 
planetary atmosphere. The entropy production due to the ir- 
reversibility of the processes occurring within the climatic 
fluid is called the material entropy production, S ma t and can 
be written in general terms as: 



mat 



-I 



F V; 



pdV 



(8) 



where the first term on the right hand side is he contribu- 
tion coming from the dissipation of kinetic energy, and the 
second term is related to the transport of heat (in sensible 
and/or latent heat forms) across gradients of the temperature 
field. One can prove that the entropy production associated 
with the dissipation of kinetic energy is the minimal value 
of the entropy production compatible with the presence of a 
Lorenz energy cycle with average intensity W, 

S m ,n= I ^ pdV (9) 
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We can define the irreversibility parameter a measuring 
the excess of entropy production with respect to the mini- 
mum, which results from the heat transport down the tem- 
perature gradient: 



Srr 



a 



Sn 



fyei/TpdV 



(10) 



which is zero if all the production of entropy is due to the 
unavoidable viscous dissipation of the mechanical energy. 
The parameter a introduced above is related to the Bejan 
number Be (Paoletti et al. 1989) as Be = a + 1. 

3 Model Setup 

Bracketing the multi- stability properties of a complex sys- 
tem by resorting to numerous runs of a numerical model 
with slightly altered values of one or more parameters is in- 
deed a computationally expensive exercise. The matter be- 
comes problematic if we are considering a system as com- 
plex as a planetary atmosphere and we are treating a two- 
dimensional parametric space. In what follows, we focus 
our attention on an Earth-like planet, and investigate the im- 
pact of changing the irradiance S of the parent star, and of 
changing the opacity of the atmosphere (modulated by the 
CO2 concentration). Our goal is to find the region in the 
(5, [CO2]) space supporting multistability for the climate 
system, and study the transitions between such states. For 
each point in the (S, [CO 2]) space where multistability is 
found, we refer to the two coexistent states as Warm (W) 
and Snowball (SB) states. Outside the region of bistability, 
the climate state will be found in either a W state or a SB 
state. 

The numerical simulations are performed 
with the Planet Simulator (PLASIM). PLASIM 
(Fraedri ch et al. 2005) is a fast-running climate 
model of intermediate complexity, freely available at 
http : / / www . mi . uni-hamburg . de/plasim. Its 
dynamical core solves the primitive equations for vorticity, 
divergence, temperature, specific humidity and the log- 
arithm of surface pressure using the spectral transform 
method (Eliasen, Machenhauer & Rasmussen 1970; Orszag 
1970) with semi-lagrangian advection. Interaction between 
atmosphere and radiation is dealt with simple but realistic 
longwave ([(Sasamori 1968)"] ) and shortwave ( [(Lacis and 
Hansen 1974)} radiation schemes. The treatment of the 
solar forcing accounts for both seasonal and the diurnal 
cycle. Unresolved processes taking place at the subgrid 
scales as moist (Kuo 1965; 1974) and dry convection, cloud 
formation and large scale precipitation (Stephens 1978; 
Stephens et al. 1982; Slingo & Slingo 1991), latent and 
sensible heat boundary layer fluxes, horizontal and vertical 
diffusion (Louis 1979; Louis et al. 1981; Laursen & Eliasen 
1989) are parameterized. The model is coupled to a 50-ra 
slab ocean which contains a thermodynamic sea-ice model. 
All simulations are performed at axial tilt angle (obliquity) 



6 = 23°. In Table [T] there is a summary of the physical and 
orbital parameters used in these simulations. 

Table 1 List of parameters and symbols. 



parameter/ symbol 


explanation 


value 


n E 


Earth's rotation rate 


7.29 • 10" b rad _i 


y 


orbital year 


360 x 8.64 x 10 4 s 


s 


obliquity 


23° 


a 


planet's radius 


6300 km 


Cd 


specific heat of dry air 


1004 J kg" 1 K" 1 




specific heat of mixed layer model 


4180 Jkg^KT 1 


9 


gravitational acceleration 


9.81 ms~ 2 


Pw 


ocean water density 


1030 kg m 3 


hml 


mixed layer depth 


50 m 


V 


pressure 




V 


velocity 




V 


volume 




p 


density 





The model is run at T21 spectral resolution (approxi- 
mately 5.6° x 5.6°) with 10 vertical levels. While this res- 
olution is relatively coarse, it is expected to be sufficient, 
for slowly rotating and phase locked planets, for obtaining 
a reasonable description of the large scale properties of the 
atmospheric circulations. A higher horizontal resolution is 
indeed needed at higher rotation rates (e.g. ft > 4£l E ) since 
the size of baroclinic structures (Rossby deformation radius, 
see, e.g. |Holton 2004| scales as We remark that pre- 
vious analyses have shown that using a spatial resolution 
approximately equivalent to T21 allows for obtaining an ac- 
curate representation of the major large scale features of the 
climate system and of its global thermodynamical proper- 
ties (Jones et al. 2005; Pascale et al. 2011). 

Finally, we consider values of S ranging from 1165 
Wm~ 2 to 1510 Wm~ 2 , and values of [CO 2] ranging from 
90 ppm to 2880 ppm. 

We wish to anticipate a caveat about the appropriateness 
of the model adopted in this study. The radiation scheme 
adopted by PLASIM has been devised to deal with Earth- 
like CO2 concentrations, and in this range it is reasonably 
accurate. While in this study we do push the radiation code 
beyond the usually considered conditions, we are still con- 
sidering CO2 concentrations within one order of magnitude 
of the present terrestrial conditions, so that we still retain 
confidence in our model. The effect of considering ultra- 
high values of CO2 concentrations has been studied in detail 
in (Pierrehumbert 2005), where it is discussed under which 
conditions using a radiative model designed for Earth cli- 
mate conditions can be problematic. 

4 Bistability, hysteresis and regime 
boundaries in the [S, C0 2 ] space 

4.1 W and SB states 

In the following we shall refer to the parametric plane 
(SXCO2 ]) as the CS space. The exploration of the CS space 
proceeds as follows for each considered value of [CO2]: 
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Fig. 1 Surface temperature (K) as a function of S and the [CO2]. The lower SB (a) and upper W (b) manifolds are 
shown. The transition SB— and W^SB are marked with dashed arrows, starting from the S s b w and S ws b tipping 
regions, respectively. 



1. the model is run to a W steady state for 100 years with 
S equal to 1510 Win 2 ; 

2. S is decreased by a small amount and the model run is 
continued until a steady state is reached; 

3. step 2 is repeated until S is reduced to 1 165 W m -2 ; the 
point of W^SB transition is noted down; 

4. the reverse operation is then performed with S increased 
step by step, up to the value of 1510 Wm -2 , each time 
allowing the system to reach a steady state; the point of 
SB^W transition is noted down. 

Using systematically this procedure, it is possible to 
identify two lines in the CS space where the SB^W and 
W— )>SB transitions occurs. We refer to these two lines as 
Ssbw ([CO 2]) and S s bw ([CO 2]). They have the property that 
for values of S larger than S s b w , only W states are possible, 
while for values of S smaller than S ws b only SB states are 
realized. Empirically, one find a rather simple parameriza- 
tion for the S s b w and S ws b lines: 

>^sbw dsbw log 10 [CO 2 ] + C ahw (11) 

Swsb = Uwsb loglo[C^2] + C ws b (12) 

where a sbw ~ a wsb ~ —70 Wm -2 , C s b w ~ 1630 Wm -2 
and C wsb ~ 1440 Wm -2 for the transition SB^W and 
W— )-SB respectively and [CO2] is expressed in ppm. We 
can estimate the width A# of the bistable region as: 

A B = C sbw - C wsb * 200Wm- 2 . (13) 

The globally averaged mean surface temperature, T s in 
the CS space is illustrated in Fig. [T] Two climatic regimes 
can be identified, corresponding to the two leaves of the 
function T 3 (S, [CO2])). The overlap of the two leaves is 



responsible for the multi stability properties, and the bor- 
der of the region where both leaves are defined correspond 
to the tipping points of the system. S s b w and S ws b are re- 
ported as dashed lines in Fig. [I] and the direction of tipping 
is also indicated. Interestingly, in the CS space, isotherms 
are parallel to the S s b w and S W8 b,so that within a very good 
degree of approximation the SB— ^W occur at T s = 242 K, 
while the W^SB take place at T s = 254 K. The bistability 
properties of the climate system are therefore "blind" to the 
mechanism of forcing and depend only on the temperature. 
One must note that surface temperature range 242 K-254 K 
is not permitted by the climate system. Across the explored 
parametric space, the temperature range on the SB and W 
manifolds are 212 K-242 K, and 254 K-326 K respectively. 
The climate sensitivity of the W state is much higher than 
that of the SB state because the SB is almost entirely dry, so 
that the positive water vapour feedback is not active and the 
surface temperature difference between the two manifolds 
ranges between 40 K and 60 K, respectively, for identical 
values of S and [CO2]. 



4.2 Transitions and parametrizations 

There is much more to say than what is portrayed in Fig. [T] 
i.e. the W and SB states differ much more than just in sur- 
face temperature. In (Bo schi et al. 201 3)] we have presented 
an extensive account of the thermodynamical properties of 
the W and SB states in the CS space. An extremely use- 
ful fact is that it is possible to establish empirical relations 
connecting thermodynamical and dynamical properties of 
a planet to a given thermodynamic quantity, which can be 
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more easily determined at experimental level, so that we can 
effectively reduce the 2D CS space to a simpler ID space. 

Specifically, one can establish approximate empirical 
laws of the form T SB (S, [C0 2 ]) « T^ B (T a (S, [C0 2 ])) 
and T w {S,[C0 2 ]) * T^(T S (S, [C<9 2 ])), where T is a 
thermodynamical property such as, e.g. entropy production, 
and the lower index refers to whether we are in the SB or 
W state. For a given quantity T, the empirical laws T^ S B and 

will in general be different. This result is quite interest- 
ing in a classical perspective of climate dynamics, where 
it is customary to parameterise large scale climate prop- 
erties as a function of the surface temperature, especially 
when constructing simple yet meaningful models ( Saltzman 
(2002)). This also suggests that, in fact, the surface tempera- 
ture, well beyond its obvious practical relevance, is, loosely 
speaking, a good climatic state variable, i.e. it contains a 
great deal of information on the physical state of the sys- 
tem. 

Similarly, it is possible to establish parametrized laws 
expressing to a good degree of approximation all the main 
thermodynamic quantities describing the non-equilibrium 
steady state of the climate system as a specific functions 
V T S E B (T E ) and Y$(T E ) of the emission temperature T E . 
This result is even more relevant in terms of both theory 
and applications. First, we find an entirely non-trivial way 
to connect non-equilibrium properties to the quantity -T E - 
that quintessential^ describes the climate system as a zero- 
dimensional, energy balance system. Second, we have a 
way to connect a quantity which can be readably observed 
also in exoplanets - T E - because only globally integrated 
energy data are needed - to quantities which are related 
to horizontal and vertical gradients of the thermodynamic 
fields of the planetary atmosphere, and thus cannot be di- 
rectly observed. 

We will discuss below these empirical parametrizations. 
Of course, one may well wonder whether the existence of 
such parametrizations is robust and, going more into detail, 
whether the empirically obtained functional relations are ro- 
bust with respect to changes in the external parameters of 
the system, like the orbital ones. Of course, one cannot per- 
form a complete study, but results are extremely encourag- 
ing at least considering a very important orbital parameter, 
the rotation rate of the planet ft. In Fig. [2] we present a scat- 
ter plot of the emission temperature T E versus the globally 
averaged surface temperature T s for ft* = 1 (magenta) and 
Q* = 0.5 (black), where ft* = ft/ft#, with ft# the current 
terrestrial rotation rate. Such a range is much larger than 
anything ever experienced by our planet, but still far able to 
encompass the case of slowly rotating planets. We will deal 
with this case in Sec.|5] 

Figure [2| shows that there is a clear monotonic relation 
between the two temperature indicators, as one could have 
guessed, and that the 242K — 254K gap in the allowed 
values of T s discussed in Fig. [T] corresponds to the gap 
233K — 238K in T E . Moreover, the (approximate) func- 
tional relation between T s and T E as well as the gaps do not 




Fig. 2 Relationship between emission temperature T E 
(K) and surface temperature T s (K) for ft* = 1 (magenta) 
and ft* = 0.5 (black). 



depend substantially on the choice of the rotation rate of the 
planet. In this figure we also introduce a simplified way for 
representing the tipping point regions. The SB^W transi- 
tions are represented by the arrow pointing from the blue 
dot to the red dot, while the W— >SB transitions are repre- 
sented by the arrow connecting the red to the blue dot. As 
we see, the projection of the CS space into a ID space is 
successful also for the critical states of the system. 

Figure [3] shows that the emission temperature is a good 
predictor of the sea-ice fraction. The SB state, defined by 
T E < 233 K, features an ocean where the sea-ice fraction is 
unity. The W state has a less trivial behavior: we can iden- 
tify a subset of states, characterized by T E > 260 K, where 
sea ice is absent, and, therefore, the ice-albedo feedbacks is 
switched off. For 238 K < T E < 260 K, there is a roughly 
linear monotonically decreasing relation between the sea- 
ice fraction and T E . It is interesting to note that transition 
from the W to SB states occurs at a critical value of the 
sea-ice fraction of about 0.5, in broad agreement with previ- 
ous results by |(Voigt et al. 20lT)| The freezing of the planet 
leads to the sea-ice fraction to unity, whereas the reserved 
SB— )-W transition has even more dramatic effects, because 
we go from a completely frozen to a virtually ice free ocean. 
This elucidates even more clearly the strength of the ice- 
albedo feedback. As we see, there is barely detectable dif- 
ference between the simulations referring to ft* = 1 and 
those referring to ft* = 0.5. 

In Figs. [4|[7] we present much less trivial empirical ther- 
modynamic relations, by expressing thermodynamic quan- 
tities as approximate functions of the emission temperature 
T E . As can be guessed by comparing with Fig. [2] exactly the 
same conclusions can be drawn by plotting the same quan- 
tities against T s , which actually provides cleaner functional 
relations ["(Boschi et al. 201 3) [ nonetheless, we choose T E as 
predictor variable because of its astrophysical relevance. 

Figure [4] shows that one can robustly predict the value 
of the material entropy production, S mat given the value of 
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Fig. 3 Sea ice fraction (in %) vs. surface temperature T s 
(K) for ft* = 1 (magenta) and ft* = 0.5 (black). 



T E . We also find that the approximate functional relation 
is weakly dependent on the rotation rate, as the clouds of 
black and magenta dots are closely entwined both in the W 
and SB states. As a reference, one may consider that in the 
present climate, the Earth's material entropy productions is 
about 50 x lQ- 3 Wm 2 K- 1 ( Pascale et al. 2011; Boschi et 
al. 2013). The values of S mat are monotonically increasing 
with Te- In the SB state, the entropy is exclusively gen- 
erated by dissipation of kinetic energy and by irreversible 
sensible heat transport, because the planet is almost entirely 
dry. As for the W manifold, the main contribution to entropy 
production comes from latent heat due to large scale and 
convective precipitation. The presence or lack of a signifi- 
cant hydrological cycle changes entirely the entropy budget 
of the planet. In the bistable region, the range of S m at is 
(10,19) 10~ 3 W m~ 2 K- 1 and (34,62) 10~ 3 W m^K" 1 
for the SB and W respectively, therefore a factor of 3 larger 
in the W regime with respect to the SB regime. Moreover, 
there is a range of values of S mat - from 19 to 34 10 _3 W 
m -2 K _1 - which is not allowed by the system. The pres- 
ence of a large gap and the very pronounced difference be- 
tween the values of S ma t m the W and SB states confirms 
that Sm a t may be a better indicator than temperature for 
discriminating the SB and W states as already discussed in 



Lucarini et al. (2010b)| and is of relevance also regarding 
the definition of habitability conditions, because large val- 
ues of Sm a t provide strong indication of the presence of an 
active hydrological cycle. 

The radical difference between the W and SB states be- 
comes more apparent when looking into measure of the 
ability of system to produce mechanical work starting from 
available potential energy. The Carnot efficiency of the sys- 
tem r] (Fig. [5} decreases abruptly with Te in the W mani- 
fold, thus implying that warmer climates are characterized 
by smaller temperature differences. The main reason for this 
is that the transport of water vapour acts as very efficient 
means for homogeneizing the temperature across the system 




210 220 230 240 250 



Fig. 4 Material entropy production S ma t 
K _1 ) vs. emission temperature T E (K) for 
genta) and ft* = 0.5 (black). 



(10- 3 W m 2 
ft* = 1 (ma- 



(Peixoto & Orrt 1992; Lucarini & Ragone 2011; Lucarini et 
al. 2011). This matches the fact that the system is charac- 
terised by very strong irreversible processes, as described 
by the very large values of entropy production realized in 
these conditions (Fig. [4]). Consequently, the value of the pa- 
rameter of irreversibility a increases as conditions becomes 
warmer and warmer (not shown); see also |Lucarini et al.| 
(2010b), When looking at the intensity of the Lorenz en- 
ergy cycle W (Fig. [6]), we discover that it is weakly depen- 
dent on Te (but monotonically decreasing with T E around 
the present climate) for ft* = 1, because the strongly re- 
duced efficiency is compensated by the stronger absorption 
due to the water vapor feedback. The ft* = 1 and ft* = 0.5 
cases do not feature a quantitative agreement comparable 
to what shown in previous indicators, especially regarding 
the intensity of the Lorenz energy cycle W, as this is di- 
rectly related to the dynamical regime of a planetary atmo- 



sphere, in turn strongly influenced by ft Pascale et al. 2013 



Nonetheless, the qualitative agreement is rather good. The 
only detectable difference is that in the W state, the value 
of W is more clearly positively correlated with T E in the 
case of slower rotation. The likely reason for the presence 
of a stiff er climate for ft* = 1 is that baroclinic instability, 
more relevant for faster rotation, is very strongly damped by 
efficient mixing of temperature. 

When considering the SB state, we observe a strength- 
ening of the dynamics of the climate system with increas- 
ing surface (or emission) temperature. This is evident when 
looking at the dependence of the strength of the Lorenz en- 
ergy cycle W (Fig. ^ and of the Carnot efficiency of the 
system r] (Fig. [5} on the emission temperature . The value 
of W reaches its largest value just at the warm boundary 
of the SB manifold. Similarly, r] increases monotonically 
with increasing temperature. This due to the fact that in the 
SB state warmer climates correspond to atmospheric condi- 
tions with reduced stratification (Bosc hTet al. 2013)] which 
allows the development of stronger large scale atmospheric 
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Fig. 5 Carnot efficiency vs. emission temperature Te (K) 
for ft* = 1 (magenta) and ft* = 0.5 (black). 




230 240 250 260 270 



Fig. 6 Lorenz energy cycle strength W (W m 2 ) vs. 
emission temperature Te (K) for ft* = 1 (magenta) and 
ft* = 0.5 (black) 



motions. The agreement is rather good for the ft* = 1 and 
ft* = 0.5 cases. 

As has been seen, the W^SB transition is associated 
to a large decrease of T s , Te, and S mat . The large drop 
in T s is related to the fact that most of the material en- 
tropy production is associated with the hydrological cycle 
that, in the quasi-dry SB states, is almost absent. Nonethe- 
less, this does not say much about the processes leading 
to the transition between the two states, or better, describ- 
ing how one of the attractors disappears. The efficiency r] 
(Fig. [5} features, instead, a more interesting behavior, with 
its maximum values just before the W^SB transition. We 
find that each transition is associated to a notable decrease 
(more than 30%) of the efficiency of the system, and the 
closer the system gets to the transition in the CS space, the 
larger is the value of the efficiency (in Fig. [5] we can see 
r] vs. T s for ft* = 1 and ft* = 0.5). This can be inter- 
preted as follows. If the system approaches a bifurcation 
point, its positive feedbacks become relatively stronger and 
the negative feedbacks, which act as re-equilibrating mech- 
anisms, become less efficient. As a result, the differential 
heating driving the climate is damped less effectively, and 
the system is further from equilibrium, since larger temper- 
ature differences are present. At the bifurcation point, the 
positive feedbacks prevail and the circulation, even if rather 
strong, is not able to cope with the destabilizing processes, 
and the transition to a state on the other manifold is realized. 
The new state is more stable and closer to thermodynamical 
equilibrium. The decreased value of the efficiency looks to 
be the marker of this property. It is interesting to note that, 
at the transitions, the values of the efficiency are practically 
equal: t]w^sb ~ Vsb,w ~ 0.03. Moreover, we observe 
that r] saturates in the very cold regime of SB states and very 
warm regime of W (Fig. [5} at rj aa t,w ~ Vsat,SB ~ 0.012. 
This leaves us with an open question: is the result on the 
decrease of the efficiency at the bifurcation points in both 
W— )>SB and SB^W transitions of more general relevance? 
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Fig. 7 Meridional heat transport index MET (in PW) vs. 
emission temperature Te (K) for ft* = 1 (magenta) and 
ft* = 0.5 (black). 

As last element of the non-equilibrium properties of the 
climate system, we discuss the relationship between Te and 
the meridional energy transport MET (Fig. [7]). We observe 
that in the SB regime there is a weakly positive relation- 
ship between Te and the transport. As discussed above, this 
is related to the fact that warmer conditions allow for a re- 
duced static stability of the atmosphere. Instead, for a vast 
range of values of Te in the W regime, the transport is al- 
most insensitive to Te (and T s ), as discussed thoroughly 
in (Caballero & Langen 2005; Donohoe & Battisti 2011; 
2012). The rigidity of the climate system has been attributed 
by Donohoe & Battisti (2011; 2012) to the fact that the 
meridional heat transport is mainly determined by planetary 
albedo and thus atmospheric composition rather that surface 
albedo and sea-ice coverage. However, such a mechanism 
ceases to exist when the ice albedo feedback becomes in- 
effective because of the disappearance of sea-ice above a 
threshold value of T E ~ 260 K. Above this value, we ob- 
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serve a steep monotonic increase of the transport with tem- 
perature, because changes in the latent heat transport are 
mainly responsible for this behavior. This agrees with the 
idea that the dynamics and sensitivity of a warm planet is, 
in some sense, dominated by the hydrological cycle. In this 
latter case, the skill of the two temperature quantities in pa- 
rameterising the MHT is comparable. Indeed, the details of 
the transfer functions will depend on various specific prop- 
erties of the planetary system under investigation. The focus 
here is on the fact that our results suggests that it is possi- 
ble to define such empirical relations. As in most previous 
cases, there is overall a very good agreement between the 
ft* = 1 and ft* = 0.5 cases; the stronger positive corre- 
lation between Te and MHT in the W state can be inter- 
preted similarly to what done when looking at W, since a 
stronger atmospheric circulation tends to have a more in- 
tense MHT. 



5 How does rotation affect bistability? 

The results presented so far have been obtained making spe- 
cific choices such as the position and size of the continents, 
the radius of the planet, the eccentricity of its orbit, the or- 
bital tilt and its rotation rate, the atmospheric composition, 
among others. Obviously, it is hard to account for such a 
vast multidimensional parametric space, and one must come 
to terms with the fact that only partial explorations are pos- 
sible (and future investigations will deal with different axes 
of such parametric space). A first parameter which is rel- 
evant in determining the dynamical and thermodynamical 
properties of the atmospheric general circulation is the rota- 
tion rate of the planet (Williams 1998a; 19988b; Hunt 1979; 
Pascale et al. 2013 ). In the previous section, we considered 
two values of the rotation rate, ft# and ft#/2. Our results 
have clarified that only minor differences exist between the 
two cases. However, the explored range falls short of explor- 
ing all possible or at least all observed planetary conditions. 
In particular, many exoplanets for which habitability condi- 
tions are debated are phase-locked as a result of tidal forces, 
so that they always face the parent star with the same hemi- 
sphere (Dvorak 2008; Kastings 2009; Seager 2010; Perry- 
man 2011), and their atmospheric circulation is driven by 
the temperature difference between the bright, warm side 
(facing the parent star), which receives radiation, has a pos- 
itive radiative budget at the top of the atmosphere, and redis- 
tributes heat towards the dark, cold side, which emits long- 
wave radiation to space ( M erlis & Schneider 2010| ). 

Such a forcing/dissipation scheme is conceptually sim- 
ilar with respect to the case of fast rotating planets like 
Earth, with the crucial difference that the boundary con- 
ditions are basically time-independent. Atmospheres and 
circulations of phase-locked planets are currently inten- 
sively investigated (|Joshi, 2003[|Heng , Frierson & Phillipps 
2011aHHeng, Frierson & Phillipps 2 011b[|Merris & Schnei-| 



very different dynamical and thermodynamical properties, 
for the basic reason that the scaling of the underlying evo- 
lution equation is entirely different from the case of fast 
rotating planets. This causes technical problems when at- 
tempting to adapt a model devised for describing terres- 
trial or Earth-like conditions for describing the circulation 
of slowly rotating planets. In our case, we have extensively 
rewritten the PLASIM's code (Fraedri ch et al. 2005)] so that, 
in a rather unusual way in the context of climate modeling, 
the evolution equations are not scaled with respect to given 
constants. This is the first step towards the construction of a 
fully flexible modeling suite for super-Earths. 

We investigate the bistability of a slowly rotating Earth- 
like planet by using a modified model set-up, where, for 
sake of generally, we choose the case of zero eccentricity, 
we remove the land- sea contrast typical of our planet, and 
an ocean-only surface is considered (so-called Aquaplanet 
set-up). We perform simulations considering the case of ex- 
tremely low rotation (ft* = 1/180, 1/120, 1/90, 1/30) plus 
the case of 1:1 phase lock condition. This corresponds to 
selecting ft* = 1/360; note that, as shown in Table [TJ the 
length of the orbital year is chosen to be equal to 360 usual 
terrestrial days. For each value of ft*, the CO2 concentra- 
tion is kept at 360 ppm, while the solar constant S is slowly 
decreased from 1500 Wm~ 2 to 1000 Wm~ 2 , and increased 
again from 1000 Wm~ 2 to 1500 Wm~ 2 . The parametric 
exploration follows the same steps as described in Section|3] 
and allows for detecting hystereses in the climatic variables, 
which can be related to the multi- stability of the system. 

Since in this case the goal is to detect potential varia- 
tions in the structural properties of the climate system re- 
sulting from changes in ft*, we focus on a single climatic 
variables. In Fig. [8] (a) -(d) we show the dependence of the 
globally averaged surface temperature T s as function of the 
solar constant S for some of the considered values of ft*. 
The red line marks the climatic states obtained when S is 
decreased 1500 Wm~ 2 — » 1000 Wm~ 2 , whereas the blue 
line describes the states obtained when S is increased 1000 
Wm~ 2 — >> 1500 Wm~ 2 . When comparing the sets of runs 
performed with ft* = 1/30 - Fig. [8(a)]- and ft* = 1/120 - 



Fig. |8(b)| -, we observe that the values of T s in the W and SB 
states are rather similar. An important difference emerges, 
though: the width of bistability A b (compare definition in 



der 2010) 7 In general, it is reasonable to expect that phase- 
locked and, more generally, slowly rotating planets feature 



Eq. 13 ) is greatly reduced for a more slowly rotating planet, 
going from ~ 150Wm~ 2 to ~ 90Wm~ 2 . The amplitude of 
the bistability interval of S is further reduced by decreasing 
ft*, up to ft* = 1/180, whic h is t he largest value of ft* for 
which Ab ~ Wm -2 (Fig. 8(c)). Correspondingly, in this 
case, near the tipping points the slope of T s as a function 
of S tends to infinity both on the left and on the side of the 
critical value S cr u ~ 1170 Wm~ 2 . 

Hence, we can introduce ft* ri£ ~ 1/180 as critical value 
for the parameter ft*, defining the occurrence of a phase 
transition in the planetary system. The bifurcation graph is 
reported in Fig. [9] where the width of the bistability region 
A# is plotted against ft*. 



www.an-journal.org 



© WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



10 



Lucarini & al.: Habitability and Multistability in Earth-like Planets 



sa->w 

W-*SB 







1200 1300 1400 1500 



(a) 



(b) 



SB->W 
• W->SB 





■ W*SB 




1200 1300 1400 1500 



(c) 



(d) 



Fig. 8 Surface temperature T s (K) vs solar constant S (fc" 2 ) for (a) Vt* = 1/30, (b) Vt* = 1/120, (c) Q* = 1/180 
and Q* = 1/360. 



For Q* < &l rit , monostability is detected and, as ex- 
pected, no divergence in the derivative of T s with respect 
to S is found, while a rather regular monotonic dependence 
dT s /dS > is found. See in Fig. 8(d)| the special case of 
1:1 phase locked planed (Q* — 1/360). In this case, steady 
states are characterised by equatorial sea-ice building up in 
the dark side of the planet and ice-free bright side; look- 
ing at the relationship between sea-ice fraction and S (not 
shown), one finds a very regular monodic dependence, with 
no threshold effects of the sort seen in Fig. [3] In the consid- 
ered range of S, the sea-ice fraction is always below unity 
and above zero. We have that the ice-albedo feedback be- 
comes weaker - to the point that bistability is lost - when, no 
matter the intensity of the incoming radiation, the equatorial 
region spends enough time on the bright side of the planet 
to allow for melting of sea ice, and, conversely, the polar 
regions spend enough time in the dark side of the planet to 
allow for freezing of surface water, unless we consider very 
low or high values of the constant S. We guess that other 
parameters of great relevance in this context are the obliq- 



uity of the planet and the relationship between the length 
of the year, which we have kept fixed to the usual one in 
these simulations, and the time scale describing the rate of 
relaxation of the temperature of the ocean's surface ( [Pascale 
|et al. 2Q13| ). 

Therefore, when no bistability is present, it is possible to 
relate one-to-one the solar constant to the globally averaged 
surface temperature (or to the emission temperature) and no 
gap in the possible values of the globally surface tempera- 
ture is found. In other words, once we lose the global in- 
stability related to the ice-albedo feedback, the system be- 
comes rather boring or, alternatively, much more easy to in- 
terpret. Habitability conditions are also more easy to assess. 



6 Conclusions 

In this review we have examined the bistability properties 
of an Earth-like planet using a set of recently developed 
thermodynamical diagnostics, which allow for defining the 
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Fig. 9 Bifurcation graph of the planetary circulation: 
width of the bistable region A B (Wm~ 2 ) vs. ft*. 

fundamental non-equilibrium properties of any planetary at- 
mosphere, ranging from material entropy production, effi- 
ciency of the climate engine powering the atmospheric cir- 
culation, and linking them to more classical dynamical in- 
dicators such as the intensity of the Lorenz energy cycle. 

We have used and modified according to our needs the 
open- source modeling suite PLASIM, an intermediate com- 
plexity global terrestrial climate model, by allowing for an 
extensive parametric exploration of rather diverse planetary 
conditions, in terms of amount of incoming stellar radiation, 
opacity of the atmosphere (modulated by the CO2 concen- 
tration), and rotation rate of the planet. 

Within a wide parametric space, which includes the 
present conditions, the climate is multistable, i.e. there are 
two coexisting attractors, one characterised by warm con- 
ditions, where the presence of sea-ice and seasonal snow 
cover is limited or altogether absent (W state), and one char- 
acterised by a completely frozen sea surface, the so-called 
snowball (SB) state. As well known, this fact has paleocli- 
matological relevance, but this is not the main direction of 
our study. 

For all considered values of [CO2] (from 90 to 2880 
ppm) the width of the bistable region is about 200 Wm -2 
in terms of the value of the solar constant, and its position 
depends linearly on the logarithm of the [CO 2] and shift- 
ing by about 15 Wm -2 per doubling of CO2 concentration. 
The W state is characterized by surface temperature being 
40 K - 60 K higher than in the SB state. In the W states, the 
material entropy production is larger by a factor of 4 (order 
of 40 - 60 lO^Wm^K" 1 vs. 10 - 15 lO^Wm^K" 1 
with respect to the corresponding SB states. The bound- 
aries of the bistable region are approximately isolines of 
the globally averaged surface temperature or of the emis- 
sion temperature, and in particular, the warm boundary, be- 
yond which the SB state cannot be realized, is characterized 
by the vanishing of the permanent sea-ice cover in the W 
regime. 

The thermodynamical and dynamical properties of the 
W and SB states are largely different. In the W states, the 
hydrological cycle dominates the dynamics and latent heat 



fluxes contribute most in redistributing the energy in the 
system and to the generation of material entropy produc- 
tion. The SB state is eminently a dry climate, with entropy 
production mostly due to sensible heat fluxes and dissipa- 
tion of kinetic energy. The response to increasing tempera- 
tures of the two states is rather different: the W states fea- 
ture a decrease of the efficiency of the climate machine, as 
enhanced latent heat transports kill energy availability by 
reducing temperature gradients, while in the SB states the 
efficiency is increased, because warmer states are associ- 
ated to lower static stability, which favors large scale atmo- 
spheric motions. The entropy production increases for both 
states, but for different reasons: the system become more ir- 
reversible and less efficient in the case of W states, while 
stronger atmospheric motions lead to stronger dissipation 
and stronger energy transports in the case of SB states. 

A general property which has been found is that, in both 
regimes, the efficiency increases for steady states getting 
closer to tipping points and dramatically drops at the tran- 
sition to the new state belonging to the other attractor. In a 
rather general thermodynamical context, this can be framed 
as follows: the efficiency gives a measure of how far from 
equilibrium the system is. The negative feedbacks tend to 
counteract the differential heating due to the stellar insola- 
tion pattern, thus leading the system closer to equilibrium. 
At the bifurcation point, the negative feedbacks are over- 
came by the positive feedbacks, so that the system makes 
a global transition to a new state, where, in turn, the nega- 
tive feedbacks are more efficient in stabilizing the system. 
On a more phenomenological note, the transition from the 
W to SB states occurs at a critical value of the sea-ice frac- 
tion of about 0.5. This agrees with previous findings (Voigt 
and Pierrehumbert 2011). After the transitions, the sea-ice 
fraction becomes unitary. When considering the reverse tip- 
ping point, the transition is even more dramatic: the sea-ice 
fraction changes abruptly from unity to virtually zero. 

We have shown that empirical functional relations 
are found between the main thermodynamical quantities 
and globally averaged surface temperature of the emis- 
sion temperature, allowing for expressing the global non- 
equilibrium thermodynamical properties of the system in 
terms of parameters which are more directly accessible. Al- 
though this method requires further investigation in order 
to delimit its range of applicability, it suggest a methodol- 
ogy to infer information about the atmospheric dynamics of 
exoplanets (which would be otherwise unaccessible), also 
because we have discovered that the transfer functions are 
rather robust with respect to large variations of orbital pa- 
rameters such as the planetary rotation rate. 

In the last part of this work we have explored the dynam- 
ical range of slow rotating and phase locked planets, where 
the length of the day and the length of the year are compa- 
rable. We have clearly found that there is critical rotation 
rate below which the multi- stability properties are lost, and 
the ice-albedo feedback responsible for the presence of SB 
and W conditions is damped. The bifurcation graph of the 
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system suggests the presence of a phase transition in the 
planetary system. As such, critical rotation rate corresponds 
roughly to the phase lock 2:1 condition. More specifically, 
the width of the bistable region is gradually reduced up to 
zero for ft* « 1/180. For ft* < 180 the two attractors 
associated with the W state (polar sea-ice caps and tropi- 
cal ice-free ocean) and the SB state (globally covered by 
sea-ice) merge and only one attractor exist, corresponding 
to a totally different climate (equatorial sea-ice in the dark 
side and ice-free ocean in the bright side). In particular, if 
an Earth-like planet were phase locked 1 : 1 with respect to 
its parent star, only one climatic state would be compatible 
with a given set of astronomical and astrophysical param- 
eters. We plan to extend the investigation of phase-locked 
planetary conditions by exploring the impact of changing 
the length of the orbital year, thus performing a thorough 
analysis of the thermodynamical properties of the the cir- 
culation patterns described by along the lines of Merlis & 
Schneider (2010). 

The results discussed in this paper support the adoption 
of new diagnostic tools based on non-equilibrium thermo- 
dynamics for analysing the fundamental properties of plan- 
etary atmospheres and pave the way for the possibility of 
practically deducing fundamental properties of planets in 
the habitable zone from relatively simple observables. Fu- 
ture investigations will analyse more systematically how ro- 
bust our findings are with respect to changes in relevant or- 
bital parameters. 
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